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Abstract. The simple supersymmetric approach recently used by Dutt, Gangopadhyaya, and Sukhatmc 
for spherical harmonics is generalized to Jacobi equation, including also the intermediate Gegenbauer case. 

Resumen. Un procedimiento supcrsimctrico simple recientemente usado por Dutt, Gangopadhyaya y 
Sukhatme para los harmonicos esfericos, es generalizado a la ecuacion de Jacobi, incluyendo tambien el caso 
intermedio de Gegenbauer. 
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I. INTRODUCTION 

During 1971-1974, supersymmetry (SUSY), a more general (graded) Lie algebraic realization of symme- 
tries, involving combinations of commutation and anti-commutation relations has been discovered as a sort of 
necessary unifying algebraic ingredient of microscopic bosonic and fermionic degrees of freedom jjj. In 1981, 
Witten |Q introduced the simplest SUSY realization within nonrelativistic quantum mechanics (SUSYQM) 
as a toy tool for studying SUSY breaking in quantum field theories. SUSYQM may be seen as a special 
property of the factorizations of the one-dimensional (ID) Hamiltonians. Namely, one can find a pair of 
Hamiltonians (usually called bosonic and fermionic partners) having the potentials connected by a single 
Riccati solution, known as the superpotential. Indeed, if one fixes to zero the ground state energy, one can 
easily factorize any ID Hamiltonian as follows 

Hl ^- c2 ^ + Vl=A ' A > (1) 

where A? = —c-^ + W(x),A = c-^ + W{x) and c = for nonrelativistic quantum mechanics. The 

W function occurring in the factorizing operators is the solution (superpotential) of the following Riccati 
equation 

Vi = W 2 - cW' . (2) 

The SUSYQM scheme, as initiated by Witten, means to construct the "fermionic" Hamiltonian obtained by 
reversing the order of the factoring operators, which reads 

H 2 = AA^ = -c 2 -^+V 2 , (3) 

where 

V 2 = W 2 + cW' . (4) 

One can show that the energy eigenvalues, the wave functions, the S matrices, the reflection and transmission 
coefficients of the supersymmetric partners Hi and H 2 are related in a simple manner through the application 
of the factorizing operators [Q . In particular, the energy spectra are identical except perhaps for the bosonic 
ground state (for the so-called unbroken SUSY case). In 1983, Gendenshtein || noticed that in many 
important problems there is a very useful relationship between the SUSY partner potentials Vi t 2, known 
as shape invariance (SI). He showed that if the superpartners are similar in shape and differ only in the 
parameters that appear in them, i.e., in mathematical terms, if one writes down a SI condition of the form 

V 2 (x;a 1 ) = V 1 (x;a 2 )+R(a 1 ) , (5) 

where ai is a set of parameters, and a 2 is a function of oi, such that the remainder R(a\) is independent of 
x, one can easily derive the energy eigenvalues and eigenfunctions of any potential fulfilling the SI condition. 
Thus, SI is a valuable parametric connection within SUSYQM. 

The coordinate transformation procedure is well known in the framework of SUSYQM (see section 5 in 
j|). Essentially, one starts with a ID Schrodinger equation and tries to obtain a better known one by per- 
forming a coordinate transformation. This may be considered as a particular case of various transformations 
of homogeneous linear differential equations of second order involving either the dependent variable, or the 
independent one, or both, which are recognized as basic methods in the mathematical literature; the reader 
is directed to the textbooks of Szego and Zwillinger B. However, using SUSYQM, one can exploit further 
SUSY symmetries, like the aforementioned SI, to get some properties of special functions in an illuminating 
way. This is precisely the case of a recent paper of Dutt, Gangopadhyaya and Sukhatme (DGS) 0, where 
they used the change of the independent variable, to recast the associated Legendre equation in the cele- 
brated Schrodinger equation with the — sech potential, which enjoys the SI property. Next, they put to work 
this concept to derive some known features of spherical harmonics in a simple way. 

In the following, after briefly repeating the DGS case of spherical harmonics, we generalize their scheme 
to the Jacobi equation, presenting also the intermediate Gegenbauer case. 
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II. SPHERICAL HARMONICS 

The equation for the associated Legendre polynomials 



de 2 



cote 



dy 

de 



1(1 + 1) 



m 
sin 2 (9 



y = Q 



(6) 



can be transformed into a Schrodinger eigenvalue equation by a mapping function e = f(z) that can be 
found from the condition of putting to nought the coefficient of the first derivative. The result is = f = 
2arctan(e 2 ). This mapping is equivalent to the replacement sin# = sechz and cos# = — tanhz. The range of 
the variable z is the full real line — oo < z < oo. We notice that 6(z) — -| +gd(z), where gd(z) is the so-called 
Gudermannian or hyperbolic amplitude function [Q. The associated Legendre equation is transformed in 
one of the best known exactly solvable, SI Schrodinger equation 



<Pv_ 

dz 2 



1(1 + l)sech 2 z 



-m 2 v 



(J) 



for which the algebraic SUSY scheme can be readily applied. Since I is always an integer for the common 
spherical harmonics, Eq. (7) is moreover a reflectionless one (i.e., the calculated reflection coefficient is zero). 
The energy eigenvalues are known to be of the type E n — —(I — n) 2 , with n = 0, 1, 2...N, where N is the 
number of bound states the —1(1 + l)sech 2 potential well can hold, which for spherical harmonics is equal 
to /. Since E n = —to 2 , one gets n = I — m. The eigenfunctions |v„(z;0)) can be obtained by repeated 
application of the creation operators A\z; k) = (— 4z + fctanhz), k s [0, m — 1], as follows 



\v n (z;0)) oc 



fe=0 



k—m—1 



\v (z;m)} , 



(3) 



where the ground state wavefunction is \vo(z;m)) = sech m z. The connection between the v functions and 
the associated Legendre polynomials is given by P ; m (tanhz) w vi^ m (z; I). 

We notice that the mapping / and the Schrodinger equation (7) were already known to Infeld and Hull 
[HI, who used them in their factorization method. 



III. JACOBI POLYNOMIALS 

A more complicated case is that of Jacobi polynomials which fulfill the differential equation 
(1 - x 2 )y" + [/3-a-(a + (3+ 2)x]y + n(n + a + (3 + l)y = . 



(9) 



By employing the notations 2(m + 1) = a + f3 + 2 , n(n + a + (3 + 1) = (n 1 — m)(n' + m + 1), (3 — a = 7. 
Eq. (9) turns into 

(1 - x 2 )y" + "fy - 2(m + l)xy + (n - m)(n + m + l)y = (10) 
and with v = (1 — x 2 ) m l 2 y one gets 



(1 - x 2 )v" - (2x - i)v' + n(n+l) 



m — jmx 



1-x 2 



v = 



or in spherical polar coordinates 
d 2 



a v / 
h 1 

de 2 \ 



cote 



7 \ dv 



smt 



)de + l n {n+1) - 



m — 7m cos f 



sm 



v = 



(11) 



(12) 



We already see that Eq. (12) presents deviations from the standard associated form. With the change of 
variable 6 = f(z) and the notations P(f) and Q(f) for the coefficients of the first derivative of v and of v, 
respectively, in Eq. (12), one can get the following form 



dz 2 



f'P(f)- y] ^ + Q(f)f' 2 v = 



(13) 
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Again the coefficient of the first derivative is put to nought, leading to the differential equation 



f P(f) 



which can be solved by the substitution u — f = leading to 

du 



df 



= uP{f) 



(14) 



(15) 



with the solution v = J P{m . Since J P{9)d6 = M ta ffi fl % ] we have £ = tan s ^"g /2) ■ Thence, z = 
J tan 7 (6'/2)csc0(i0 + const. We have found that a convenient choice for the constant of integration is const = 
——. Using g — tan(0/2), one gets z = J g 1_1 dg — — . Thus, 

C tan^e/2) _ 1 , Q 

Z = < 7 7 

I ln[tan(0/2)] 7 =0. 



The angular variable can be written as 



_ / 2arctan[( 7 z + I) 1 / 7 ] 7 ^ 
2arctan(e z ) 7 =0. 



One can show that the mapping / for 7 7^ is equivalent to the replacements sin 9 = ^ z+1 ^-i/^ + ^ z+1 -ji/-y 

sechw and cos 6 = j^fg^ffegff 
get an equation of the form 



tanhw, where w = M7£+il_ p or the general case 7 7^ one can 



d^ 



n (n + 1) 



sech 2 ( ln(7Z+1) ) m 2 + 7 m tanh( ln(7Z+1) ) 



(72 +1) 5 



(jz + I) 2 



v = 



(16) 



It can be interpreted as a Schrodinger equation at zero energy. One set of solutions v are of the type 

V n = [1 - tanh 2 ( ln (7*+l) °+£ pa.fl ( tanh ( ln( 7 »+l) ) ) > where pa,/3 arg j acobi po l ynomia l s . 

Up to now we did not use the a, asymmetry in the function v. If we do that by using v = (1 — x)"/ 2 (1 + 
x) lS / 2 y in Eq. (10), one gets 



(1 — x 2 )v — 2xv + n (n + 1) — 



5 + ex 



1-x 2 



v = , 



where <5 = m(2m + 1) + 7 2 and e = —7(3™ + 1). In spherical coordinates (x = cos9) Eq. (17) reads 



d 2 v dv 

-dF 2+coW Te 



n (n +1) 



5 + e cos ( 



sin 2 



u = , 



(17) 



(18) 



which though not in the standard associated Legendre form has the advantage that the change of variable 
9 — 2arctan(e z ) can be used to get 



dz 2 



[n (n + l)sech 2 z + etanhz]i> = Sv . 



A simplified form containing only the tanh function is the following 

d 2 v 



dz 2 



[Atanh 2 z + 2£tanhz]v = Cv 



(19) 



(20) 



where A — —n'(n + 1), B = — ^ ( 3 " +1 ) ; C = —n(n + 1) + m(2m + 1) + 7 2 , and it is supposed that 
B < n 2 B. As discussed by DGS and by other authors [^0| such a Schrodinger equation corresponds to 
the motion of an electron in a Coulomb field in the presence of an Aharonov-Bohm potential and a Dirac 
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monopole potential. It is also a Rosen-Morse II problem which is known to be shape invariant Q. The 
energy eigenvalues and eigenfunctions are in our notations as follows || 



E n = n (n + 1) - (n - nf - — - , (21) 

[n — n) z 

v n = (VI - tanhz)"'-" +r (\/l + tanhz) n '""- r i 3 ™'- n+r '' l '- n - r (tanhz) , (22) 

where r = , B . 

n —n 

IV. GEGENBAUER POLYNOMIALS 

The Gegenbauer polynomials (x) are polynomial solutions of the ultraspherical equation 

(1 - x 2 )y" - (2q + l)xy + pip + 2q)y = . (23) 

The ultraspherical case can be considered as the particular 7 = Jacobi case in the formulas 9-16 of the 
previous section, and at the same time as a slightly generalized spherical harmonics case of the first section. 
One gets immediately the shape invariant, exactly solvable Schrodinger equation 

v = -m' 2 v , (24) 

where n = p + q — \ and m = q — 5. For a full equivalence to Eqs. (6) and (7) q should be half integer. 
The energy eigenvalues of the potential V(z) = —n'(n' + l)sech 2 .z are E n = —{p! — n) 2 , (n = 0, 1, 2, ...N), 
where N is the number of bound states in the potential well, and is equal to the largest integer contained in 
n'. The eigenfunctions v n (z;n') are obtained by applying the factorization (creation) operators A^{z;ai) — 
(— + at tanh z), where ai — n' — i, onto the ground state wave function vo(z;a n ) = sech™ ~ n z = sech 9 ~ 5 z, 
i.e., 

v n (z- n') w A\z- n')A*<{z; n' - l)A^(z; n' - 2)...A t (z; n'-n+ l)sech"'- Tl z , (25) 

which is an analog of Eq. (8). In the case of this section n = n' — m' = p, and therefore we are dealing with 
associated Legendre functions of the type v n = P™ (tanh 2;). The original Gegenbauer polynomials can be 

written formally as C^(x) oc (1 — tanh 2 z)^ £ ~P < ^ + ^_ 1 (tanh z); the proportionality is used since we are not 
addressing the normalization issue here. Although this relationship appears to be known, see jnj, we found 
it in a simple way. The n'(n' + l)sech 2 z potential is reflectionless if and only if n' is an integer. Since p 
should be an integer, the reflectionless property requires again a half integer q. 

Gegenbauer polynomials are connecting higher-dimensional spherical harmonics to the usual two- 
dimensional ones. In the theory of hyperspherical harmonics they play a role which is analogous to the 
role played by Legendre polynomials in the usual three-dimensional space [^2|. Therefore, we believe our 
discussion to be useful in such a context as well. 
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